ON THE NUMBER OF FACETS OF THREE-DIMENSIONAL 
DIRICHLET STEREOHEDRA III: FULL CUBIC GROUPS 



PILAR SABAPJEGO AND FRANCISCO SANTOS 



Abstract. We are interested in the maximum possible number of facets that 
Dirichlet stereohedra for three-dimensional crystallographic groups can have. 
The problem for non-cubic groups was studied in two previous paper by 
D. Bochis and the second author. This paper deals with "full" cubic groups, 
while "quarter" cubic groups are left for a subsequent paper. Here, "full" 
and "quarter" refers to the recent classification of three-dimensional crystallo- 
graphic groups by Conway, Delgado-Friedrichs, Huson and Thurston. 

Our main result in this paper is that Dirichlet stereohedra for any of the 
27 full groups cannot have more than 25 facets. We also find stereohedra with 
17 facets for one of these groups. 



1. Introduction 

Tiling problems are among the most classical problems in discrete geometry, yet 
they include many basic ones that are still open. See [H Chapter 4] for a recent 
survey. A fundamental question, explicitly asked in [U p. 177], [HI p. 959], and 
[2"U] . and related to Hilbert's 18th problem [IB] is: 

What is the maximum number of facets of a convex polytope P 
such that K 3 can be tiled by congruent copies of PI 

The fact is that there is no known upper bound for the number of facets of such 
a polytope. Only in the special case where the symmetry group of the tiling acts 
transitively over the tiles, a global upper bound is known. In this case, the tiling 
is called stereohedral and the tile is a stereohedron: 

Theorem 1.1 (Dclonc A stereohedron in Euclidean d- space with respect to 
a symmetry group G cannot have more than 2 d (a + 1) — 2 facets, where a is the 
number of aspects of G. 

The group G C Isom(R d ) that acts transitively on a stereohedral tiling is nec- 
essarily a crystallographic group. That is, it contains a full-dimensional discrete 
lattice T = 1 d as a subgroup of finite index. The quotient group G/T is usually 
called the point group of G, and its cardinality is the number of aspects of G. Equiv- 
alently, the number of aspects of G equals the number of translational lattices in 
which a generic orbit of G decomposes. In dimension three, crystallographic groups 
can have up to 48 aspects. Hence, the above "fundamental theorem of the theory 
of stereohedra" implies that three-dimensional stereohedra cannot have more than 
390 facets. 

On the other hand, the stereohedron with the maximum number of facets known 
so far has "only" 38 of them. It was found by P. Engel in 1980 (see [5] and [141, 
p. 964]), using a computer search, and it tiles space under the action of the cubic 
group I4i32, with 24 aspects. 
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It is also interesting to observe that "unbounded stereohedra" , that is, un- 
bounded polyhedra that tile M. d by congruent copies under the action of a dis- 
crete transformation group, can have arbitrarily many facets, even in dimension 3. 
See [11], or Lemma 4.2 in [2]. 

A special case of stereohedral tiling is the Voronoi diagram of an orbit of points in 
a crystallographic group. The stereohedra obtained are called Dirichlet stereohedra 
or plesiohedra. The fact that Engcl's 38 face stcrcohedron is actually a Dirichlet 
stereohedron seems to indicate that Dirichlet stereohedra can have as many facets 
as general stereohedra (although there is no proof of this). Since Dirichlet stere- 
ohedra are much easier to analyze, in previous papers the second author together 
with D. Bochis, has initiated an exhaustive study of the number of facets of Dirich- 
let stereohedra for the different three-dimensional crystallographic groups. They 
proved: 

• Within the 100 crystallographic groups which contain reflection planes, the 
exact maximum number of facets of Dirichlet stereohedra is 18. Within 
cubic groups containing reflection planes the maximum number of facets is 
eight pQ. 

• Within the 97 non-cubic crystallographic groups without reflection planes, 
no Dirichlet stcrcohedron can have more than 80 facets [2] . The same paper 
shows an explicit construction of Dirichlet stereohedra with 32 facets for a 
tetragonal group. 

For cubic groups without reflection planes (there are 22 of them), Bochis, and 
Santos were only able to prove an upper bound of 162 facets Here and in [TH] 
we improve this bound, and hence the general upper bound for the number of 
facets of three-dimensional Dirichlet stereohedra. For technical reasons, we do this 
separately for "full groups" (this paper) and "quarter groups" (paper [Hj). Among 
other things, the proof we present here for full groups is "computer-free", while the 
proof for quarter groups is not. 

Our work takes great advantage of the new classification of crystallographic 
groups developed by Conway et at. in [5]. To enumerate cubic groups, they use 
the fact that every cubic group contains rotations of order three in the four diag- 
onal directions of the orthonormal lattice Z 3 . They call odd subgroup of a cubic 
crystallographic group G the subgroup generated by all these rotations and show 
that: 

Theorem 1.2 (Conway et al. [5]). For every cubic crystallographic group G: 

(1) Its odd subgroup is one of two possible groups R\ or i?2- 

(2) Ri and R2 are normal in Isom(M 3 ), hence G lies between Ri and its nor- 
malizer Af(Ri), for the appropriate i. 

In particular, the classification of cubic groups is equal to the classification of 
(conjugacy classes of) subgroups of the finite groups N(Ri)/R\ and A/"(i?2)/^2, 
which are of orders 16 and 8. (The latter is isomorphic to the dihedral group D%; 
the former contains D$ as an index two subgroup) . The odd subgroups R\ and R2 
differ in that R\ contains rotation axes of order three that intersect one another 
while all rotation axes of order three in R2 are mutually disjoint. Conway et al. call 
full groups those with odd subgroup equal to R\ and quarter groups those with odd 
subgroup equal to i?2- (The reason is that R2 is a subgroup of R\ of index four. 
That is, quarter groups contain only one quarter of the possible rotation axes). 
There are 27 full groups (14 of them without reflection planes) and eight quarter 
groups (none with reflection planes, since A/"(i?2) does not contain reflections). 

The main result in this paper is: 
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Table 1. Our bound for the number of facets of Dirichlet stere- 
ohedra of full cubic groups without reflections. The six groups 
marked with an asterisk are symmorphic. Their Dirichlet stcre- 
ohcdra were completely described by Engel (1981), and have at 
most 17 facets. 

Theorem 1.3. Dirichlet stereohedra for full cubic groups cannot have more than 
25 facets. 

More precisely, for each of the 14 full groups without reflections we obtain the 
upper bound shown in Table [1] Here and in the rest of the paper we use the 
International Crystallographic Notation for crystallographic groups, see [T7j. Let 
us remind the reader that for cubic groups with reflections, Theorem 2.4 in |T] gives 
a bound of eight facets. 

Although, as far as we know, our work represents the first significant improve- 
ment to Delone's global upper bound on the number of facets of general Dirichlet 
stereohedra, there are at least the following previous classifications of special Dirich- 
let stereohedra: 

• Engel [9] completely classifies Dirichlet stereohedra for the symmorphic 
cubic groups. G is called symmorphic if it is the scmidirect product of its 
translation subgroup and its point group. Equivalently, if there is an orbit 
of G on which the translational subgroup acts transitively. 

Symmorphic cubic groups are obviously all full. There are nine of them 
with reflections six without reflections. Engel found that their Dirichlet 
stereohedra have up to 17 facets. 

• Koch [TS] completely classified Dirichlet stereohedra for cubic orbits with 
less than three degrees of freedom. She found 117 types of stereohedra, 
with up to 23 facets. 

• Dress et al. [7] completely classified the stereohedral tilings whose symm- 
metry group acts transitively on the two-dimensional faces. There are 88 
types of tilings, and 7 types of tiles. 

Remark 1.4. In the previous papers in this series we always assumed that the 
base point p for the orbit has trivial stabilizer for the crystallographic group G 
under study. We believed this was not a loss of generality in the sense that every 
crystallographic orbit with non-trivial stabilizer for a group G is also an orbit with 
smaller stabilizer for another group G' . That is actually not true, but we prove in 
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the appendix that there are only four groups G in which it fails, all of them cubic 
and without reflections. These cases were previously found by Fischer [12 H3] by 
an exhaustive study of cubic orbits with less than three degrees of freedom. 

We also show in the appendix that in the full case these special orbits produce 
stereohedra with at most 12 facets. In the quarter case we refer to the classification 
of Koch mentioned above, which gives a bound of 23 facets. 

2. The structure of full cubic groups 

The odd group R\ of full cubic groups is a crystallographic group of type F2Z (in 
International Crystallographic Notation). It is generated by rotations of order three 
in all the lines with directions (±1, ±1, ±1) and passing through integer points. In 
particular, the points where two or more axes meet form a body centered lattice 

7 = {(0,0,0),(i i ±)} + Z. 

To understand full groups, the Delaunay tesselation of this lattice is specially help- 
ful. Each Delaunay tetrahedron has two opposite dihedral angles equal to 7r/2, and 




Figure 1. The body-centered cubic lattice /. The Delaunay tetra- 
hedra fall into two classes, according as they can be moved to co- 
incide with the tetrahedron on the left or right, by a motion that 
preserves the labeling of the points in I. 

the other four dihedral angles equal to 7r/3. The latter lie along rotation axes of 
order three and the former on rotation axes of order two, parallel to coordinate 
axes. Since every edge has an even number of incident Delaunay tetrahedra, the 
Delaunay tesselation is 2-colorable. The group R\ acts transitively and with trivial 
stabilizer on the Delaunay tetrahedra of a given color. 

Another way to define the coloring is as follows. The action of R\ divides the 
lattice / into four orbits, namely: 

F, := {(a, b, c) e I : a + b + c - | £ 2Z}, i = 1, 2, 3, 4. 

Each orbit is a face centered cubic lattice. F2 and F4 together form a primitive 
lattice, and F\ and F3 form another. 
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Each Delaunay tetrahedron has a vertex in one of each of the Fj's, and two De- 
launay tetrahedra have the same color if and only if the isometry that sends one to 
the other and keeps the labels of vertices is orientation-preserving. In Figure Q] lat- 
tice points are labelled according to their four orbits, and two Delaunay tetrahedra, 
which belong to different "colors" . have been shaded. 

We fix, for the rest of the paper, the following "base Delaunay tetrahedron" T: 

T := conv({(l/2, -1/2, 1/2), (1, 0, 0), (1/2, 1/2, 1/2), (0, 0, 0)}. 

T is one of the two shaded tetrahedra in Figure [1] We call its four vertices «i , v-i , 
V3 and i>4, in the order we have given them, so that v% G Fi. We also introduce the 
following notation: 

• Let Ti, i = 1, ... ,4 be the tetrahedron adjacent to T and sharing with it 
the triangle opposite to w,. Let v[ the vertex of Tj not in T . 

• Similarly, let , with i j, denote the neighbor of T, and sharing with it 
the triangle opposite to Vj. Its vertices of are Vk, vi, v[ and a certain 
vertex v'/j, where {i,j, k, 1} — {1, 2, 3, 4}. 

Observe that: 

• If i — j is odd, then T, TJy and Tji are related by the rotation of order three 
on the axis passing through vu and vi. 

• If i — j is even, then v'/j = v[. That is, T\z = T31 and T42 = T24. Moreover, 
Tij is related to T by an order two rotation on the axis passing through Vk 
and vi . The rotation of order four on this same axis relates the two of them 
with Ti and Tj. 

The normalizer Af(Ri) of Ri is the symmetry group of the lattice /. It acts 
transitively over all the Delaunay tetrahedra and with a stabilizer of order eight, 
isomorphic to the dihedral group Dg (this stabilizer is the symmetry group of the 
Delaunay tetrahedron itself). Hence, the index of R\ in Af(Ri), that is, the order 
of the quotient group Af(Ri)/Ri, is 16. 

Knowing Af(Ri) and R\ we can give the whole list of all full groups. They 
are pictured in Table [5J We classify them with two parameters: the order s of the 
stabilizer of T, and a number m = or m = 1 depending on whether the group acts 
on the Delaunay tesselation of the lattice I with one or two orbits of tetrahedra. 

A group G with m = is completely characterized by the stabilizer of the base 
tetrahedron T. Indeed, G is generated by i?i together with this stabilizer. Hence, 
such groups are in bijection to the eight (conjugacy classes of) subgroups of Dg, 
the symmetry group of T. The lattice of them is drawn in Figure [H For each 
group we have drawn a projection of T in the vertical direction and, inside it, the 
orbit of a generic point. The dots representing an orbit are drawn some black and 
some white, depending on whether the corresponding point lies above or below 
the horizontal plane at height 1/4, which bisects T. The corresponding full cubic 
groups are listed in the first column of Table with the same picture. The three 
groups that contain reflection planes are shaded. 

For a group G with m = 1, we first look at the transformations in G that send 
Delaunay tetrahedra to Delaunay tetrahedra of the same i?i-orbit. They form a 
subgroup H of index two in G, and H is one of the eight groups of the previous 
paragraph. Then, G is generated by H together with any single transformation 
p G G that sends T to one of its neighbors, say T4. Let Ht denote the stabilizer of 
T in H , and let px be p composed with the mirror reflection on the plane through 
Vi, V2 and U3 (that is, the unique transformation that sends T4 to T, keeping the 
labels of vertices). Let Gt be the group generated by Ht and pr- Then, G is 
completely characterized by Ht and Gt and the conditions on Ht and Gt are 
that H T <G T < D s and \G T ■ H T \ equals 1 or 2. 
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Figure 2 . The subgroups of symmetries of the base tetrahedron 
T. 

To picture these groups, in Table we show the image in T4 of the orbit of 
H pictured in T. The rule of thumb is that the two orbits, when superimposed 
by composing one with the mirror symmetry on the plane containing v\ , 1*2 and W3 
should together form another orbit of a subgroup of Dg. This happens automatically 
when s = 8 or s = 4. For s = 1, it is equivalent to the transformation between the 
two orbits being an involution of T. For s = 2, it is necessary but not sufficient 
that the transformation be an involution. 

Depending on the particular H, there may be one or more possibilities for G. 
There is always at least one, in which p is the mirror symmetry itself (and the 
involution is the identity). That is the first group shown after H in each row of 
Table [21 and it is shaded since it obviously contains reflections. The other ones 
contain reflections if and only if H contains reflections. 

3. A FIRST UPPER BOUND 

Let now G be one of the 27 full groups, let T be a Delaunay tetrahedron of the 
lattice I, and p a point in T, whose Dirichlct stereohedron we want to study. We 
assume that p has trivial stabilizer for the action of G on Gp. That is to say, that 
the only motion in G that fixes p is the identity. In the appendix we treat the case 
of non-trivial stabilizers, and correct statements made about this in the previous 
papers [U [2] . The main reason for this assumption is the following lemma: 

Lemma 3.1. If p has trivial stabilizer, then a small perturbation of it cannot 
decrease the number of facets of the Dirichlet stereohedron Votq p (p) . That is, the 
maximum number of facets for a Dirichlet stereohedron with trivial stabilizer is 
achieved when p is generic. 

As a consequence, from now on we assume that p is generic, in particular that 
it lies in the interior of T. 
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Table 2. The 27 full cubic groups. 



Proof. Let pi G Gp is a neighbor of p in the Voronoi diagram. We claim that when p 
is perturbed to a position p' sufficiently close to p the point p\ £ Gp' corresponding 
to p\ will be a neighbor of p' . 

For this, let x G R 3 be a point whose distance to p and p\ is equal, and smaller 
than the distance to any other point in Gp. By continuity, a small perturbation of 
p will make p' and p' x still be the two orbit points closest to x, except perhaps not 
equidistant to it. Assume, without loss of generality, that after the perturbation 
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we have d{x,p') > d(x,Pi). Then, moving a; in a straight-line towards p' makes 
d(x,p') decrease faster than the distance to any other point in Gp' . Hence, when 
x reaches a position equidistant to p' and p\ , the distance to all other orbit points 
will be bigger, which proves that p' and p\ are neighbors in V otq p i . □ 

Lemma 3.2. The Dirichlet stereohedron Voyg p (p) is fully contained in the (non- 
convex) polyhedron obtained as the union of T and its four neighbors T\, T2, I3 
and T4. Moreover, no facet ofVoTG p (p) lies in the boundary of this polyhedron. 

Proof. We start with the following simple observation: if a group G contains a 
rotation p of angle a and axis on the line I, then Vor g p (p) is fully contained on 
the dihedral sector of angle a with edge I and centered around p. Indeed, this 
dihedral sector equals Vor s(p) where S C Gp is the orbit of p under the cyclic 
group generated by p. 

In our case, since the six edges of T are rotation axes for every full group, 
VoiG P (p) will be contained in the intersection of the six corresponding dihedral 
sectors. For each rotation of order three on edges of T, say the one on the edge 
v±V2, the sector of angle 2tt/3 in question is in turn contained in the half-space 
delimited by the plane containing the vertices vx, V2, v' 3 and v\. The intersection 
of the four half-spaces is an infinite square prism circumscribed to T. Similarly, 
the rotations of order two on the axes and V2V±, cut this infinite square prism 
by two sectors of 37r/2 with dihedral angle at the edges v\v 3 and V2V4. The non- 
convex polyhedron obtained in this way coincides with the union of T and its four 
neighbors T\, T2, T3 and T4. 

For the "moreover" observe that p being interior to T implies that the half-spaces 
and dihedral sectors can be taken open, except for the corresponding axes. □ 

Lemma 3.3. Let q 6 Gp be a point whose Dirichlet stereohedron VorQ p (q) shares a 
facet with Vor Gp(p)- Then, q lies in the interior of one of the following 15 Delaunay 
tetrahedra of I: T , its four neighbors Ti, or the ten additional neighbors of the latter 
T 

Proof. By definition, p and q arc neighbors if and only if Voro P (p) and Vor<3p(<z) 
share a facet. Lemma l3~2l implies that this facet (except perhaps some of its edges) 
is contained in the interior of the union of T and its neighbors. The same lemma 
applied to q, implies that q lies cither in T, or in a neighbor of T, or in a neighbor 
of a neighbor of T. □ 

Since each Delaunay tetrahedron contains at most eight points of Gp (because 
Gp C j\f(Ri)(p), which has exactly eight points in each tetrahedron), this lemma 
gives already a global upper bound of 8 x 15 — 1 = 119 facets for Dirichlet stere- 
ohedra of full cubic groups. Even more, full groups that contain eight points per 
tetrahedron necessarily have reflection planes. So, every full group without reflec- 
tion planes has at most 4 orbit points in each tetrahedron, which gives a bound of 
4x15 — 1 = 59. But a more refined bound can be given using the following ideas: 

• First, we will give bounds for each group G separately, using as parameters 
the numbers s(G) and m(G) introduced in the previous section. We do this 
to identify which groups are worth of a closer study. For example, in case 
m(G) = then the list of 15 tetrahedra given in Lemma \'6. 3 1 goes down to 
only 11: T and the ten "neighbors of neighbors" Ty of T. 

• Second, the following lemma allows us to do almost as if the number of 
"neighbors of neighbors" of T was only six and not ten. 

Lemma 3.4. Let G G Isom(K. 3 ) be a crystallographic group and p G G a rotation. 
Let p G K 3 and let q a point of the orbit of p, Gp. Let S be the orbit of q by 
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the rotation p. Then, at most two points of S are neighbors of p in the Voronoi 
diagram of Gp, namely, those which make the smallest dihedral angle with p, in 
both directions, when seen from the axis of p. 

Proof. Let S — S q be the orbit of q by the rotation p, and let S p denote the orbit 
of p. Clearly, (S p U S q ) C Gp implies Votg p (j>) Q Vors p uS q {p) an d, in particular, 
points of S q that are neighbors of p in Votq p must also be neighbors in Vovs p us q - 

We now claim that if a point q' £ S q is a neighbor of p in Vors us then the 
relative interiors of the regions Vor s (p) and Vors q (q') intersect (this is actually a 
variation of Lemma 1.2 in [2] and can be stated more generally). The reason is that 
if a point x lies in the common facet of the Voronoi regions of p and q' in Vor,g us, 
then p and q' are the unique closest points to x in S p and S q , respectively. 

With this claim the proof is easy to finish. The Voronoi diagrams of both S p 
and Sq consist of dihedra with an edge in the axis of p and with angle 2tt divided 
by the order of p. The dihedron Vor<? (p) can only overlap with (at most) two such 
dihedra of the diagram Vors 9 . 

Perhaps the special case where p and q lie in the same orbit modulo (p) needs 
an extra word. In this case, S — S p — S q and the Voronoi diagram Vors pU s,, (p) 
equals the cycle of dihedra around the axis of p mentioned above. It is then clear 
that in it every cell (in particular, the one of p) has exactly two neighbors. □ 

In the following statement, we say that a full cubic group G "mixes colors" if it 
contains transformations that send black tetrahedra to white tetrahedra, and vice 
versa, in the 2-coloring of the Dclaunay tessclation of the lattice I described in 
Section H 

Corollary 3.5. Let G be a full group. Let s{G) be the order of the stabilizer of T 
in G, and let m(G) be equal to 1 if G mixes colors of Delaunay tetrahedra and if 
it does not. Then, the number of facets of any Dirichlet stereohedron for the group 
G is bounded above by 

(7 + 4m(G))s(G)+3. 

Proof. The neighbors of p in the Voronoi diagram of Gp fall in one of the three 
following cases: 

• Those in T, which are at most s(G) — 1. 

• Those in one of the four neighbor tetrahedra, which are at most 4s (G), but 
are present only if m(G) = 1. That is, they are 4s(G)m(G). 

• Those in one of the ten "neighbors of neighbors" of T. In principle the 
bound would be 10s(G), but it gets reduced to 6s(G) +4 using Lemma I5~4l 
Indeed, let S = {pi, . . . ,p s (G)} be the set of points of Gp inside the base 
tetrahedron T, with p = p\. Similarly, let Sij be the set of points of Gp 
inside the tetrahedron . Then, Sij is obtained from S by the rotation of 
order two or three that sends T to T^. In particular, Lemma [3~4l says that 
for (i,j) equal to one of (1,2), (1,4), (3,2) and (3,4), S^ and Sji contain, 
in total, at most s(G) + 1 neighbors of p in the Dirichlet tessclation of Gp: 
the two obtained from p\ by the rotation of order three, and only one of the 
two obtained from each of the other points in S. That is, the total number 
of neighbors in the ten T^-'s is bounded above by 2s(G) + 4(s(G) + 1) = 
6s(G) +4. 

Adding up, the number of neighbors is at most 

s(G) - 1 + 4s(GMG) + 6s(G) + 4 = 4s(G)m(G) + 7s(G) + 3. 



□ 
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Table [3] shows the bound given in this statement for each of the 27 full groups. 
Those that contain reflections appear in parentheses. As a conclusion, the global 
upper bound derived from this statement is 47. 



s(G) 


m(G) 


Bound 


Groups 


1 





10 
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14 
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2 
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4 
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(H^ (^3), (P£3£), (/43m) 
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Table 3. The upper bound derived from Corollary 13. 51 



4. Groups that contain rotations of order four 

Corollary 13.51 is ultimately based on the fact that every full group contains the 
rotations of order three or two on all the edges of the Delaunay tetrahedra (see the 
proof of Lemma I3~2")) . But, of course, the groups for which the bound obtained is 
worse are precisely those that have more transformations. Here and in the next 
section we take advantage of this fact by using additional rotations present in some 
of the full groups. 

In this section we look at the groups which contain the rotations of order four 
in the two edges of T with dihedral angle of tt/2. We call pis and P24 these two 
rotations (py is meant to have axis on the edge ViVj). Certainly, every group 
containing them has m = 1, since both rotations send T to one of its neighbors. 
Moreover, to locate them in Table [U we just need to observe that P13 and P24 act 
on the four labels 1, 2, 3 and 4 of points in the lattice I as the transpositions 2^4 
and 1 <-» 3, respectively. So, a group from Table contains these two rotations if 
and only if the corresponding pictures contains two specific dots, namely those in 
the picure for the group P432 (said differently, P432 is the group generated by Pi 
and the rotations pis and p24)- 

That is, the groups we are interested in are 

4-2 4-2 4-2 
P432,P432,P-3-,(P— 3— ) and (J— 3—). 

n n mm mm 

As usual, those with reflection planes are in parentheses. 

Our goal is to prove the analogs of Lemmas 13.21 and 13.31 in this case. That is to 
say, we want to specify a region guaranteed to contain Votg p (p) and derive from 
it a region guaranteed to contain all neighbors of p. To do this we introduce the 
following setting, which will be used again in the next section: 

(1) We divide the base tetrahedron T in cig ht smaller tetrahedra T A ,T B ,T C , 
T D , T E , T F , T G , and T H , via the two mirror reflection planes of T and 
their two bisectors. These eight tetrahedra are congruent to one another. 
The labels A to H are given cyclically, in the way shown in Figure [31 That 
is, T A is incident to vertex V\ and edge Uii>2, T B is incident to vertex V2 
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and edge V1V2, and so on. The same decomposition will be assumed in all 



Figure 3. The decomposition of a Delaunay tetrahedron into 
eight "fundamental subdomains" . 

the other Delaunay tetrahedra, with similar notation. For example, T^ 2 
denotes the image of T by the unique transformation that sends T to 
T12 preserving the labels of vertices or, equivalently, the sub-tetrahedron 
in T12 incident to its vertex v'{ 2 and to the edge v'{ 2 V3. Figure [3] not only 
shows the decomposition of T, but also those of its neighbors Ti, T 2 , T3 
and T4. The Tj are drawn smaller and using the same projection we used 
for T. Observe that this projection makes the vertices v[ and Vi +2 coincide 
(indices assumed modulo four). 

We call the eight tetrahedra in which each Delaunay tetrahedron is di- 
vided fundamental subdomain. This is justified by the fact that they are 
fundamental domains of the normalizer Af(Ri). In particular, they are 
permuted transitively by the normalizer Af(G) of any full group. 
(2) Using the same arguments as in Lemma 13.21 we will construct a region 
VorExt(T" 4 ) with the following properties: it is a union of several fun- 
damental subdomains and for every p in the interior of T , Votg p (p) is 
contained in VorExt(T j4 ) and with no facet of Yovq p {p) in the boundary of 



We call this region an extended Voronoi region of T A . Lemma 13.21 savs 
that the union of T and its four neighbors (which is the union of 40 funda- 
mental subdomains) would be a valid VorExt(T' 4 ). But, of course, we want 
to find one as small as possible, since this will produce better bounds. 

Lemma 4.1. Let G be a full cubic group. Let S be a fundamental subdo- 
main, obtained from T A by a certain transformation g £ Af(Ri). Then, for 
every q in the interior of S we have 





cly ¥,' 



V 4 




^ V l 



VorExt(T A ). 



Vor G<? (<?) C 5 (VorExt(T A )). 
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Proof. Since g £ Af(Ri) C Af(G), we have that Vor Gq (q) = g(Vor Gp (p)), 
where p = g^ 1 (q). This, together with Votg p (p) C VorExt(T A ), finishes 
the proof. □ 

In the conditions of this lemma, we call sr(VorExt(T A )) the extended 
Voronoi region of the subdomain S, and denote it VorExt(S'). 
(3) We call influence region of T A the union of all the fundamental subdomains 
S such that VorExt(S') and VorExt(T A ) overlap. We denote it Infl(T^). 
We have that: 

Corollary 4.2. Let p be a point in the interior of T A and let q £ Gp be a 

point whose Dirichlet stereohedron Vor Gp(l) shares a facet with Votg p (p). 
Then, q lies in Infl(T" 4 ). 

Proof. Same as the proof of Lemma 13.31 □ 

In particular, Lemma 13.31 implies that Infl(T" 4 ) will be contained in the 
union of T, its four neighbors, and the ten neighbors of the latter. 

Let us mention that this method of using fundamental subdomains, extended 
Voronoi regions, and influence regions was already used in [1] and [2], with the 
same names, although applied to planar crystallographic groups. 

Figured] shows the Extended Voronoi region obtained for the groups that contain 
the two rotations of order four. Our method to compute it is to start with the 

3 





Figure 4. The extended Voronoi region for groups that contain 
rotations of order four. 



union of T and its four neighbors and cut some fundamental subdomains from it 
via the additional rotations present in these groups. For example, the rotation P13 
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allows us to exclude the whole of T2 and half of T4 and the rotation P24 allows 
us to exclude the whole of T\ and half of T3. Finally, the order two rotation po 
in the axis perpendicular to both v\Vj, and V2W4, which is an element in all these 
groups, allows us to exclude the fundamental subdomains T D , T E and T F . So, the 
Extended Voronoi region consists of the following eleven fundamental subdomains: 

VorExt(T A ) = {T A , T B , T c , T G , T H , Tf, T 3 S , T 3 G , T-f , T 4 A , T 4 B , T 4 C , if}. 

The influence region is shown in Figure [5l Its computation is a bit more com- 
plicated. In principle, one should compute the 120 extended Voronoi regions for 
the fundamental subdomains contained in the 15 tctrahedra of Lemma 13.31 and 
check which ones overlap VorExt(T j4 ). But things can be slightly simplified. On 
the one hand, we have observed (a posteriori) that the influence region does not 
increase when adding T D , T E , and T F to VorExt(T A ). That is, we can act as 
if the extended Voronoi region of any given fundamental subdomain S consisted 
of the whole Dclaunay tetrahedron containing S plus half of two of its neighbors, 
namely those closest to S. With this simplification, the reader should be able to 
check quickly that the influence region consists of the 48 fundamental subdomains 
of Figure O which are three full Delaunay tetrahedra plus six halves of them. 





»13 



Figure 5. The influence region for groups that contain rotations 
of order four. 

Corollary 4.3. Dirichlet sterohedra for groups of types P432, 7432 and -P^3^ 
cannot have more than 11, 22 and 23 facets, respectively. 
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Proof. We start with the case of G = P432. If p is a point in T A , then Gp contains 
exactly one pair of antipodal points in T, and in every Dclaunay tetrahedron P; or 
Tij. Since Infl(T j4 ) is made of 12 "half-tetrahedra" , it contains exactly 12 points of 
Gp, that is, eleven plus p itself. 

Similarly, since 7432 and P^3- contain two pairs of opposite points in each of 
P, T l or T i:i , Infl(T A ) contains exactly 24 points of Gp, that is 23 plus p. In 7432 
we subtract one, because in this case we have counted points in both P 3 ^ and Pj^, 
of which only one can be a neighbor of p by Lemma 13.41 □ 

5. Groups that contain the transversal rotation of order two 

We call "transversal rotation of order two" the symmetry po of T that exchanges 
Vi <-> i>3 and V2 «-> i>4- That is, the rotation of order two on the axis perpendicular to 
the edges V1V3 and V2V4, and crossing both of them. (It is one of the transformations 
used in the previous section to cut out the extended Voronoi region). 

Equivalently, in this section we are interested in the full groups that contain 
P23 as a subgroup. The list of them, excluding the ones already considered in the 
previous section and the ones that contain reflections, is: 

P23, 123, P-3, P4 2 32, P43n. 
n 

Figure [HI shows the extended Voronoi region. The rotations used to cut it from 
the union of T and its four neighbors P, are: 

• The transversal rotation po itself, which excludes the following fundamental 
subdomains: 

fromP: T D ,T E ,T F , 

from Pi : Pf 4 , Pf , Pf , Pf , Pf , pf , 

fromP 2 : T? ,T 2 D ,T 2 E ,T 2 F . 

• The rotation of order two on the axis V1V3, which excludes the following 
fundamental subdomains from P 2 : 

rj^A rj^B rj-iC rj~>D 

• The rotation of order two on the axis V2V4, which excludes the following 
fundamental subdomains from T\ : 

rr\A rriB rr>G rr\H 

1 1 > 1 1 t 1 ! i 1 ! ■ 

• The rotation of order two on the axis perpendicular to V1V3 and V2v'^; that 
is, the conjugate of po with axis crossing P4 and P 2 . It excludes the following 
fundamental subdomains from P4: 

rpE rpF 

J-4 ,1 4 ■ 

• The rotation of order two on the axis perpendicular to viv f 3 and V2V4; that 
is, the conjugate of po with axis crossing T\ and P3. It excludes the following 
fundamental subdomains from P3: 

rpD rpE rpF 

1 3 ' 3 J 1 3 ■ 

We leave it to the reader to verify the above assertions. There are nineteen 
fundamental subdomains not excluded by any of the above, namely those 
not shadowed in Figure [H] 

We also leave out the details of the computation of the influence region. The 
interested reader, should verify that sending VorExt(P^) to be based at each of 
the 31 shadowed fundamental subdomains of Figure [7] produces a region VorExt(S') 
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Figure 6. The extended Voronoi region for groups that contain 
the transversal rotation of order two. 



that does not intersect VorExt(T" 4 ) itself. Hence, the influence region of T A consists 
of the 89 not shadowed fundamental subdomains. 

These 89 subdomains are divided as follows, according to the action of the group 
R\. Representatives for the sixteen orbits of fundamental subdomains via this ac- 
tion are, for example, the eight subdomains contained in T and the eight contained 
in any specific neighbor Tj. The numbers of them arc given in the following tables: 



Orbit 


Subdomains 


Bound on 




in Infl(T A ) 


nbr. of facets 




9 


8 


rpB 


9 


6 


rpC 


6 


4 


rpD 


6 


4 


rpE 


11 


7 




8 


7 


rpG 


6 


4 




6 


4 


Total 


61 


44 



Orbit Subdomains 
in Infl(T A ) 



Total 28 



In the case of the subdomains in T, the table contains two numbers. The first num- 
ber is the number of fundamental subdomains of that orbit contained in lnH(T A ). 
That is, a mere counting of the "white triangles" of each type in Figure [7] The 
second number is the amount by which the subdomains of a given type contribute 
to the upper bound on the number of facets. It is obtained from the first number 
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Figure 7. The influence region for groups that contain the 
transversal rotation of order two. 



by subtracting 1 to T A because of the base subdomain p, plus one to the corre- 
sponding class by the presence of the following pairs of subdomains in InflfT" 4 ). Of 
each pair, only one subdomain may produce a facet, by Lemma 13.41 



(T&,T 3 B 2 ) 

(T 2 E 3 ,T 3 E 2 ) 

{Tg,Tg) 
(T 2 H 3,T£) 



{Tg,Tg) 
(T-R,Tg) 
(TH,Tg) 

(Ti%,Tg) 
(TH,Tg) 



(T 2 F 3 ,T 3 F 2 ) 



In order to obtain the specific bound for each of the five groups we are interested 
in, we just need to add, for each of them, the bounds corresponding to the funda- 
mental subdomains that contain a point of Gp, for each specific G. The result is in 
the following table. 
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Group G 


subdomains containing 

r\Amrc r~iT ( T n 
JJUlIlLb KJL \JTJJ 


Bound on the 

IlLlIllULjl *J1 IdLLLo 


P23 
123 
P-^3 

n 

P4 2 32 
P43n 


rpA rpE rpC T~>G 
rpA rjiE rpB rp F 

2^ 

J^G 


8 + 7 = 15 
8 + 7 + 3 + 3 = 21 
8 + 7 + 4 + 4 = 23 
8 + 7 + 6 + 7 = 28 
8 + 7 + 4 + 4 = 23 



6. The group P4 2 32 

In this section we analyze in a bit more detail the group on which the methods 
so far give the worst upper bound, namely P4 2 32. Summing up, what we have said 
in the previous sections is that for every choice of a base point p in the fundamental 
subdomain T^: 

• Only the 37 orbit points lying in the following fundamental subdomains 
can possibly be neighbors of p in the Voronoi diagram of the orbit Gp: 



rpA 

J 13i 


rpB 

J 13' 


rpE 
rpE 

J 12; 


rpF 

J 13' 


rpA 

1 14i 


rpB 

J 24) 


rpE 

J 24) 

rpE 

J 21> 


rpF 

J 24) 

rpF 

J 21) 


rpA 

J 14: 


rpB 

1 14' 


rpE 

J 14> 


rpF 

± 14' 


rpA 

J 4U 


rpB 

1 41> 


rpE 

■Ml) 


rpF 

-Ml) 


rpA 

J 23i 


rpB 

23 ' 


rpE 

+23' 


rpF 

+23; 


rpA 

J 32) 


rpB 

j 32j 


rpE 

J 32) 


rpF 

J 32) 


rpA 

J 34' 


rpB 

J 34> 


rpE 

J 34i 




rpA 

J 43' 


rpB 

J 43' 


rpE 

J 43> 


rpF 

M3) 



• The count gets down to 28 since the two orbit points in T?j- and T^, with 
a; G {B,E, F} and |i — j'| + 1 2. cannot be neighbors for the same choice of 
p fLemma l3.4p . 

In what follows we investigate in further detail which fundamental subdomains 
produce a neighbor for all choices of p, or for none, and also which further depen- 
dences are there between fundamental subdomains that prevent them to produce 
neighbors simultaneously Our main goal is to lower the upper bound a bit for this 
particular group, but also this analysis gives an idea of how far our bound is from 
being tight. 

We start with the following relatively easy observations: 

Lemma 6.1. For any base point p in general position in the region : 

(1) The orbit points in T^ , T^, T43, Tq , T^, and form a Delaunay cell, 
with the combinatorial type of an octahedron in which the opposite point to 

is either T^ A or . 

(2) Similarly, the orbit points in , T^, T^, Tq ' , T^ 3 , and T^ 2 form a Delau- 
nay cell, with the combinatorial type of an octahedron in which the opposite 
point to is either or T^ 2 ■ 

(3) The four orbit points T^ , Tq , Tq , and T Q F also form a Delaunay tetrahe- 
dron. 

(4) Let 7T be the plane that bisects Tq along the mid-points of the four edges 
viVi + i. Let 7r + and ir~ be the half-spaces containing V\ (and 1)3) and V2 
(and V4), respectively. If p is in n + then Tq , Tq , T^ and form a 
Delaunay tetraedron. Otherwise, Tq- , Tq , T^ and T^ form one. That 
is, for every base point either both T-fl and T^ or both T^ and T^\ are 
neighbors. This is not exclusive, in principle the four of them could be 
neighbors simutaneously. 

Proof. In part one, the six points in question are equidistant to the mid-point of 
the segment UiV 2 , and closer to this mid-point than any other orbit point. Actually, 
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the points in T^, T^ 4 , form an equilateral triangle on the side of v\ and those 
in Tq , T^ 4 , and T 4?i form another, on the side of V2- Depending on the direction of 
the plane V1V2P, the points of the second triangle closer to p will be Tq and T® 4 or 
T fl andTj. 

The same arguments, on the edge i)\Vi, prove part 2. 

Part 3 holds since these four points are equidistant and closest to the centroid 
of T. For part (4), the four points in question are equidistant and closest to the 
midpoint of V1V3 or to the mid-point of V2U4, depending on whether p lies in the 
half-space ir + or tt~ . □ 

Corollary 6.2. For every base point in general position, the Dirichlet stereohedron 
formed by the crystallographic group P4232 has at least 11 facets. More precisely: 

• T A> T 43> T o> T 2Z, T 32> T o > and T o are always neighbors o/T A . 

• Exactly one of T^ 4 or and one of T 23 or T^ 2 is a neighbor of T§ . 

• At least one of or and at least one of T 24 or T 24 is a neighbor of 

rpA 

1 o ■ 

We have done experiments with 150 base points in general position in and 
the result is that the number of facets has always been between 14 and 17. This 
shows that the lower bound of the corollary is fairly good, while the upper bound 
of 28 can perhaps be lowered. In what follows we actually lower it to 25. 

Our experimental results are shown grafically in Figures [8j [9] and 1101 In each 
picture, the 37 fundamental subdomains that can possibly produce neighbors is 
drawn. Subdomains are colored white, grey or black according to whether the 
observed behaviour is that they always, sometimes, or never produce a neighbor. 
We do this for an arbitrary base point (Figure [HJ and also separately for base points 
chosen in the lower and upper halfspaccs 7r~ and tt + of part 4 of Lemma 16. II 

The first observation is that the statement in Corollary [62] is tight in the sense 
that all the subdomains that are white in the pictures are mentioned in that state- 
ment. In fact, the only whites that does not explicitly follows from Lemma 16.11 
is that, apparently, when p is in the lower half-space automatically and T 23 
become the opposite vertex to in the two octahedra referred to in parts 1 and 
2 of the Lemma. 

But more important for our purposes is to understand the black subdomains in 
the figures. What we experimentally observe is: 

(1) Both in the upper and lower case, only five of the eight possible subdomains 
in the tetrahedra T13 and X24 produce neighbors. More precisely, if p is in 
the lower half-space then T 24 , T 24 , and are never neighbors, while if p 
is in the upper half-space then T^, T^, and T 24 are never neighbors. 

(2) The following six subdomains never produce a neighbor: T% 4 , T43, Tf 2 , T 2l , 

j 14j j 23- 

(3) There are another four and three, respectively, subdomains, that never 
produce a neighbor for a base point in the lower and upper half-spaces. 

In our next and last result we prove part 1 of the above list, which decreases our 
upper bound for the group P4232 from 28 to 25. Proving part 2 would decrease it 

to 23 (subtract one for each pair T^ 4 T 43 , Tf 2 T 2l ) and proving part 3 would 

further decrease it to 21. 

To lower the bound we consider the subgroup of motions of G that send the 
vertical axis of To (the line passing through the midpoints of V1V3 and ^2^4) to 
itself. This subgroup can be generated, for example, by the three rotations of order 
two that stabilize To plus those on the edges v\v^ and ^2^4 of To (two of these five 
rotations are redundant, but not any two). 
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Figure 8. Experiments with the group P4232. Arbitrary base 
point. 

Let G denote this subgroup, and let Gp denote the orbit of an arbitrarily chosen 
point in the fundamental subdomain . 

In order to simplify our arguments, we slightly translate the system of coor- 
dinates that we have been using so far to one in which the vertical axis of To is 
actually a coordinate axis and with the origin in the centroid of To. More precisely, 
we take as coordinate system the one in which 

«i = (0, -1/2, 1/4), v 2 = (1/2, 0, -1/4), v 3 = (0, 1/2, 1/4), v 4 = (-1/2, 0, -1/4). 

The generating rotations of G then are on the following axes: 

• The vertical coordinate axis {x — y = 0}. 

• The diagonal horizontal axes at height zero {z = 0,x + y = 0} and {z = 
0,x-y = Q} 

• The horizontal axes in the coordinate directions at height 1/4 and/or height 
-1/4. 

More generally, G contains rotations of order four in the "diagonal horizontal" 
lines at height every integer or half-integer, and the "coordinate horizontal" axes 
at every height n/4 for an odd integer n. It also contains the two screw rotations 
of order four with axis {x = y = 0} and translation vector (0, 0, ±1/2) 
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Figure 9. Experiments with the group P4232. Base point in 
lower half-space. 



We can also relax a bit our conditions on the base point p, so that for the purposes 
of the following theorem we can consider To to be the whole strip delimited by the 
planes z = ±1/4, and divide it into eight sectors such as the following one that we 
call T A . 

T A := {(x,y,z) G K 3 : -1/4 < z < 1/4, x > 0,x + y < 0}. 

We also use the rest of labels Tq or T?j in this extended sense. 

One nice feature of this subgroup G", highlighted by our choice of coordinates, 
is that the reflection that interchanges the two half-spaces 7r + and 7r _ (now the 
reflection on the plane z — 0) is in its normalizcr. In particular, the choices of 
base point in the "upper" and "lower' part of produce equivalent Dirichlct 
tesselations. In what follows we assume that the base point is chosen in the upper 
part of T^, but what we say is true for the lower part too except there is the 
following involution that has to be applied to the names of subdomains when passing 
from an upper to a lower base point, or vice- versa: those in Tq remain the same, 
while those in T13 and T24 are exchanged as follows: 

rpA rpE rpB rpF rpE rpA rpF rpB 

J 13^ J 24: j 13^ j 24j ^13^^24: J 13 < ^- £ 24- 
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Figure 10. Experiments with the group P4232. Base point in 
upper half-space. 

Of course, only the subdomains of types A, B, E or F will contain orbit points 
when the base point is in Tq . Figure [TT1 schematically shows a vertical projection 
of these subdomains together with (part of) an orbit for this group. The figure is 
designed as an aid in the proof of the following result: 

Theorem 6.3. All Voronoi diagrams for a generic orbit of the group G are equiv- 
alent, and their Voronoi cells have nine facets. More precisely, the neighbors of 
a generic base point in the upper half of T A are the orbit points in the following 
eleven fundamental subdomains: 

rpE rpF rpB rpE rpF rpA rpE rpF rpA rpA J rpF 

^ 13> - 1 13> J i J ' J ' 24 ' 24 ' 24 ' -\- ! — ' aaa 1 + i 

where Tf = + (0, 0, 1), T A = T A - (0, 0, 1) and Tf = T F + (0, 0, 1) denote the 
translates of T A and Tq one unit up or down. 

Observe that we have three "extra" neighbors, such as , in this statement 
which do not appear anywhere else in this paper. The reason is that this is a 
neighbor in the Dirichlet tesselation of the group G' (not a crystallographic one) 
that is never a neighbor in a crystallographic full cubic group. Observe also that 
the notation is a bit misleading, since this fundamental subdomain should be 
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Figure 11. (Part of) an orbit of the subgroup G' of the group 
G = P4 2 32. 

considered "ot type E" rather than A. Indeed, the screw rotation or order two with 
translation vector (0,0, 1) is an element in the odd subgroup, and gives Tf as the 
transformed of Tf , while the translation (0, 0, 1) itself is not in the odd subgroup. 

Proof. Let p = p = (a,—(3,h) be our base point, with < h < 1/4, < a < 
/3. We introduce the following names for the other relevant points in the orbit. 
The following list also shows the fundamental subdomain they lie in. As in the 
statement, the subdomain T± denotes the subdomain T§ translated one unit up or 
down: 



P-4 ' 


= (-a,(3,h-l) 


e T E , 


9-4 = 


= (a,-f3,h-l) 




P-S : 


= (a, [3, -ft -1/2) 


G T A 


9-3 = 


= (-a,-/?, -ft - 1/2) 




P-2 : 


= 09, a, ft -1/2) 


fc 1 13' 


9-2 = 


= (-/?, -a, ft -1/2) 


GTf 3 , 


P-l z 


= (A -a, -ft) 


f= T B 
fc J J 


9-1 = 


= (-/3,a,-h) 


er F 


Po = 


(a, -(3, ft) 




9o = 


[-a, j3, ft) 


F T B 
t ^0 ' 


Pi = 


(-a, -/3,-h+ 1/2) 


e T&, 


9i = 


(a, A -ft + 1/2) 


eT 2 i, 


P2 = 


{- (3, -a,h+ 1/2) 




92 = 


(/3, a, ft + 1/2) 




P3 = 


(-/3, a, -ft + 1) 




93 = 


(/?, -a, -ft + 1) 


erf, 


PA = 


(-a,0,h + l) 


erf, 


94 = 


(a, -/3,ft + l) 


erf, 



That the neighbors of po must be contained in this list follows from the fact that 
the translations (0, 0, ±1) lie in G' . 

The structure of the orbit is a double helix: The points pi are cyclically ordered 
on the edges of an infinite prism over a (non-regular) octagon, and each point g, 
is "opposite" and at the same height as the corresponding pi . The pairs of points 
have their heights ordered by i, since < ft < 1/4. Figure [TT] shows (the projection 
of) the spiral of points pi, with p<± superimposed over p_4, indicating for each point 
the subdomain in which it lies. 
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From this double helix structure we can derive some properties of the Dirichlet 
tesselation or, rather, its dual Delaunay triangulation: 

• For every i, the tetrahedron {pi-\, qi—x,pi, qi} is a Delaunay tetrahedron. 
Indeed, its four points are equidistant and closest to the point (0, 0, i/4). 

• For every i, the tetrahedron {qi-2, qi~i,Pi,Pi+i} is a Delaunay tetrahedron. 
Here we work out the details for i = 0, the other cases being symmetric to 
it. 

Consider the empty sphere with center at the point o = (— c, — c, 0), for 
a positive c. The orbit points closest to it for c sufficiently large are clearly 
Pi and q~2- On the other hand, when c approaches zero, the closest points 
to o are those with minimum (in absolut value) height and lying in the 
half-space x + y < 0; that is, po and q_%. Moreover, it is easy to show 
that at any intermediate c the closest points to o must be one of these two 
pairs. Hence, there must be a c where these four points are equidistant and 
closest to o, hence they form a Delaunay tetrahedron. 

• Similarly, for every i, the tetrahedron {pi_2,Pi-i, Qii Qi+i} is a Delaunay 
tetrahedron. 

These three classes of tetrahedra do not cover completely the convex hull of the 
orbit (which is an infinite prism over an octagon) but almost: indeed, for each i 
these three tetrahedra form a "layer" separate the prism into two parts, so that 
every other Delaunay tetrahedron must lie in between two layers. See Figure fT2l 
where (the vertical projection of) the layers with i = and i = 1 are drawn. 




Figure 12. Two consecutive layers (three tetrahedra each) in 
VorG'p - In between the two layers lie the four tetrahedra 
q-2q-iPiP2, q-iqoPoPi, P-iPoqoqi, and p_ 2 P-i<7i<72- 

Moreover, the i — 1th and the ith layers have three edges in common, so that 
the space between them consists of four non-overlapping and triangle-disjoint tetra- 
hedra. These four tetrahedra must then also be Delaunay (and add no Delaunay 
edges except for Piqi±4 ), which finishes the description of the Delaunay tesselation. 

Hence, the list of neighbors of po in the Delaunay tesselation is 

q-4,q-3,q-2,q-i,qo,qi,q2,q3,<l4,P-i,Pi, 
as claimed. □ 

Wc remark that a "single helix" analogue of this last statement was proved as 
Lemma 4.2 in [5J Lemma 2.8] and, independently, in [10| Lemma 2.1]. 

Appendix and erratum. Orbits with non-trivial stabilizer 

In this and the previous papers in this series we have always assumed that 
the base point p for the crystallographic orbit Gp under consideration has trivial 
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stabilizer in G. That is, there is no non-trivial element g £ G with gp = p. Wc 
believed this was no loss of generality "since otherwise Gp is also an orbit of a 
proper subgroup of G" (cf. [2J p. 92]). This sentence is actually false. In this 
appendix we show that: 

• The sentence only fails in tetragonal and cubic groups, and only for orbits 
with stabilizer of order two. 

• In the tetragonal case, the sentence is almost true: every tetragonal crys- 
tallographic orbit Gp is also an orbit with trivial stabilizer for a different 
group G', except G' may not be a subgroup of G. 

• In the cubic case, there are exactly four (one-parameter families of) cubic 
orbits that are not orbits with trivial stabilizer for any crystallographic 
group. Two arise in full groups and produce stereohedra with at most 12 
facets, and two in quarter groups. For the latter, results of Koch [T5] give 
a bound of 23 facets. 

A.l. Non-trivial stabilizers have at most order two. Let I? be a stereohedron 
for a crystallographic group G (not necessarily a Dirichlet one). Suppose that the 
stabilizer of D in G is a non-trivial subgroup Ho < G but that, still, D is not a 
stereohedron for any proper subgroup of G. That is, there is no proper subgroup 
F such that {fD : / £ F} = {gD : g £ G}. 

Observe that Hq can naturally be considered a subgroup of the point group Go 
of G: every element of Hq fixes the centroid of D, which we take as the origin. 

Lemma A.l. There is no proper subgroup Fq < Go such that FqHq = Gq. 

Here F H$ = {hf : h £ H , f £ F } is not necessarily a subgroup of Go- It is a 
union of right cosets of Fo and also a union of left cosets of Hq. It is guaranteed to 
be a subgroup only if one (or both) of Fq and Hq are normal. 

Proof. Let H and F denote the preimages of Hq and Fo under the canonical quotient 
map G — > Go- Suppose that FqHq = Gq and let gD be a copy of D obtained by 
some element g £ G. We want to show that gD = fD for some / £ F. 

Let go £ Go be the point transformation corresponding to g. Since FqHq = Gq, 
there are elements ho £ Hq and fo £ Fo such that /o/io = So- Since Hq < G, we can 
consider ho as an element of G, not only a point transformation. Let / £ F be any 
element that projects to fo. Then, fhoD = fD is a translated copy of gD. That 
is, gD = tfD for some translation t. But / £ F implies tf £ F, which finishes the 
proof. □ 

Corollary A. 2. Hq is contained in all the subgroups of Gq of index two. 

Proof. If Fo < Go has index two and Hq is not contained in Fo then FoHq = Go. □ 

Theorem A. 3. If Hq is not trivial then it has order two. More precisely, either 

(1) G is tetragonal, and Hq is generated by a vertical rotation of order two; or 

(2) G is cubic, and Hq is generated by a rotation of order two in one of the 
coordinate directions. 

Proof. (1) Let P 2l , with i £ {1, 2,3} denote the point group Put differently, 

Pq, Pi and P2 are the symmetry groups of a hexagonal prism, square prism, and 
parallelepiped, respectively. The last one is contained in the other two. 

Every non-cubic point group is contained in one of these three groups. More 
precisely: triclinic, monoclinic and orthorhombic point groups are contained in P2; 
tetragonal point groups are those contained in F4 but not on P2; and trigonal and 
hexagonal are those contained in Pq but not on P 2 . So, let i be such that Go < P-n- 
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For each index two subgroup N of P-n, N n Go is either Go or an index two 
subgroup of Go- Hence, by the previous corollary, Hq is contained in the intersection 
of the following three subgroups of P-n-, of index two: The group of orientation- 
preserving symmetries, the group of transformations sending the positive vertical 
direction to itself, and the group Pi of "symmetries of the i-gonal prism" (Pi is the 
order four group generated by the mirror reflections with respect to a vertical and 
the horizontal plane). 

In the three cases this intersection is the group of order i generated by a rotation 
around the vertical axis. (Of course, the "rotation of order one" is the identity). If 
i = l then the intersection is trivial, and the case i — 2 is the one allowed in the 
statement. In the case i = 3, if Hq is not trivial (hence equal to the G3 generated 
by a vertical rotation of order three) we can apply Lemma lA.ll to the index three 
subgroup Fo := P2 R Go of Go- The three elements of C3 belong to the three 
different cosets of P2 in Pq, hence to the three of Fo m Go- 

(2) If G is cubic, then Go is contained in the symmetry group of the regular 
cube. By the corollary, Hq must be contained in the subgroups of orientation 
preserving transformations and of symmetries of the inscribed regular tetrahedron. 
The intersection of them is the group 23 of order twelve consisting of the identity, 
the three rotations of order two in coordinate directions, and the eight rotations of 
order three in the four diagonal directions. But Hq cannot contain a rotation p of 
order 3. If it did, let Fo be the subgroup of Go consisting of transformations that 
send the z axis to itself. Since 1, p and p 2 lie in the three cosets of Fj in Go, we 
have that FqHq = Go- 

So, Hq must be contained in the group 222 of order four generated by the ro- 
tations of order two in the three coordinate directions. It cannot equal that group 
either, or otherwise we could take Fj the subgroup of Go that sends the diagonal 
axis {x = y = z} to itself. Since the four diagonal axes are permuted transitively 
by the group 222, FqHq would equal Go- □ 

A. 2. Non-trivial stabilizers exist only in cubic groups. Here we assume that 
p is a base point with stabilizer generated by a rotation of order two in the vertical 
direction of a tetragonal group G. We prove that Gp is also an orbit with trivial 
stabilizer for another group G'. 

Since G is tetragonal and Gp lies in order two rotation axes, the horizontal 
projection of Gp is a square lattice. For the sake of concreteness, let us introduce 
coordinates in which this lattice is the integer two-dimensional lattice Z 2 . 

Let Gh be the subgroup of G generated by the vertical rotations of order two on 
axes that contain orbit points. Let Th < Gh be the subgroup of translations in Gh- 
Th is the group of horizontal translations by vectors with even coordinates. Put 
differently, Th splits the rotation axes in Gh into four orbits, with representatives 
in the vertical lines ?(o,o)i ^(1,0) > ^(1,1) an d ^(0,1) ■ ihx,y) denotes the vertical line 
through {x, y, 0)). 

Also, let T v be the minimal vertical translation under which the orbit Gp is 
invariant, and let c be its length. On each vertical axis, Gp splits into either one or 
two orbits of T v and it thus is an orbit of T v together with perhaps a reflection. We 
choose the vertical coordinate so that the orbit along the line ^(o,o) is symmetric 
under the change of sign in this coordinate. 

Only the three possibilities are compatible with G being a tetragonal group: 

• On all the l( XjV -)'s, the orbit points lie at the same height. Then, we let G' 
be generated by the translations with integer horizontal coordinates and 
vertical coordinates multiple of c, together with a mirror reflection in the 
horizontal plane at height zero, if needed. 
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• On all the l( x ,y) ,a the orbit points lie at different heights. This can happen 
only if G contains a screw rotation p of type 4x (or 43) on the vertical axis 
through the point (1/2, 1/2,0). The translation part of p is of length c/4. 
We let G' be generated by Th, T v , p and a rotation of order two on the 
horizontal axis passing through (0,0,0) and (1, 1,0). By our choice of the 
third coordinate this rotation keeps the orbit Gp invariant. Hence Gp is 
still an orbit of G' and it has trivial stabilizer in it. 

• On i(o,o) an( i hx,i) ^ ne orbit points are at the same height, but different to 
that of l(i.o) an d ^(0,1) (and the latter are equal to one another). If this 
happens, then either along each vertical axis there is only one translational 
orbit, or there are two but then the orbit on the axes I (0,0) an d '(1,1) is 
displaced a vertical length of exactly c/2 with respect of the orbit on in 0) 
and £(0,1) (if none of this happens, then there is no transformation of 3- 
space that sends Gp to Gp and switches the two translational orbits along 
the axis i(o,o))- I n the first case we let G' be generated by Th, T v , the 
translation with vector (0, 1, 1), and a rotation of order two in a horizontal 
axis through the points (0, 1/2, a) and (1, 1/2, a). Here, a is half the vertical 
displacement between the orbits in i(o,o) an d i(i,o)- I n the second case, we 
let G' be generated by Th, T v , and the translations with vectors (1,1,0) 
and (0,1, c/2). 

A. 3. Dirichlet stereohedra with non-trivial stabilizer, for cubic groups. 
Fischer [T2] completely classified the cubic orbits with less than three degrees of 
freedom. This includes the case of all cubic orbits with non-trivial stabilizer. His 
findings in this respect are as follows (see also [13]): 

Proposition A. 4 (Fischer). There exist exactly five cubic groups with orbits that 
are not orbits with trivial stabilizer for non-cubic groups: the quarter groups /43d 
and I^3§, and the full groups F§3 ; F4i32, and F±3%. 

Here we recover this result and show that, in the full case, Dirichlet stereohedra 
produced by these orbits have at most 12 facets. The group F^3^ is superfluous 
for us, since the "bad" orbits in it are also orbits of its proper subgroup F4i32. 

Theorem A. 5. Let p be a base point with non-trivial stabilizer in a full cubic 
group G, and suppose Gp is not an orbit with smaller stabilizer for any other crys- 
tallographic group. Then, G is one of the groups F4i32 or F|3 and the Dirichlet 
stereohedra of this orbit have at most twelve facets. 

Proof. By Theorem IA.3[ p lies in an order-two coordinate-parallel rotation axis. 
There are two types of coordinate-parallel rotations of order two in full cubic groups. 
Those on the long sides of the fundamental subdomains and those in the common 
perpendicullar line to the two long sides of one fundamental subdomain. 

If p lies in a rotation axis of the second type, then Gp is an orbit for some proper 
subgroup of G: If there is one orbit point in each fundamental subdomain, then 
Gp is an orbit with trivial stabilizer for either F23, F432 or F^3. If there are two 
(symmetrically placed) orbit points in each fundamental subdomain, then Gp is an 
orbit for either F4i32, or F%3^. See Table 2. 

1 ' an 

If p lies in a long side of the fundamental subdomain, suppose that it lies on 
the edge V1V3 of T and closer to v\ than to U3. We say that p lies "near v\ \ The 
whole orbit Gp decomposes into at most four orbits of F23, one "near v" for each 
i = 1,2, 3, 4. There are four cases, depending on the action of G on the four symbols 
1, 2, 3 and 4 of the lattice points. That is, on how G mixes the four sublattices 
of type F that make the I lattice whose Delone tetrahedra we call fundamental 
subdomains: 
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• If G does not mix the sublattices, then the orbit consists of six points near 
each center of the lattice F, lying in the two directions of the three coordi- 
nate axes from that center. This is also an orbit for the semi-direct product 
of the translational lattice F and the order six point group 3 generated by 
the inversion and a diagonal rotation of order three. 

• If G mixes the four sublattices, then the same is true except now for the 
lattice I. 

• If G mixes 1 only with 3 and 2 only with 4, then the same holds for the 
lattice P. 

• So, the only problematic case is when G mixes the sublattice 1 with one of 
2 and 4 and the sublattice 3 with the other. Looking at Table 2, we see that 
the only full groups that do this are F^3, F4i32 and F-^3%. The latter 
is discarded since then the stabilizer of p contains not only the rotation on 
the axis U1U3, but also a mirror reflection. 

In this case we consider Gp as the union of two orbits of the index 
two subgroup i*"23. The Dirichlet stercohedra for each .F23-orbit are the 
octahedra obtained glueing the four fundamental subdomain that share a 
long edge, the edge being V1V3 for one orbit and V2V4 for the other. Then, 
p has at most eight neighbors in its own i ? 23-orbit (corresponding to the 
eight facets of the octahedron), and at most four in the other F23-orbit 
(because each octahedron from one orbit intersects four from the other). 

□ 

For quarter groups, a quick look at the list of them (Figure 2 in that paper) 
easily implies: 

Lemma A. 6. If p has non-trivial stabilizer in a quarter group G and Gp is not an 
orbit for a proper subgroup of G, then G is either I43d or 7^3^. 

Proof. Of the other six quarter groups, three do not have rotations of order two 
and the other three contain the first three as subgroups of index two. □ 
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